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Picard’s Existence and Uniqueness Theorem 

Denise Gutermuth 


These notes on the proof of Picard’s Theorem follow the text Fundamentals of Differential 
Equations and Boundary Value Problerns, 3rd édition, by Nagle, Saff, and Snider, Chapter 
13, Sections 1 and 2. The intent is to make it easier to understand the proof by supplementing 
the présentation in the text with details that are not made explicit there. By no means is 
anything here claimed to be original work. 

One of the most important theorems in Ordinary Differential Equations is Picard’s 
Existence and Uniqueness Theorem for first-order ordinary differential équations. Why is 
Picard’s Theorem so important? One reason is it can be generalized to establish existence 
and uniqueness results for higher-order ordinary differential équations and for Systems of 
differential équations. Another is that it is a good introduction to the broad class of existence 
and uniqueness theorems that are based on fixed points. 


Picard’s Existence and Uniqueness Theorem 

Consider the Initial Value Problem (IVP) 

y' = f(x,y), y(x 0 ) = y 0 . 

Suppose f(x,y) and ^(x,y) are continuous functions in some open rectangle R = {(x,y) : 
a < x < b,c < y < d} that contains the point (. xo,yo ). Then the IVP has a unique solution 
in some closed interval I = [x 0 — h,x 0 + h], where h > 0. Moreover, the Picard itération 
defined by 

X 

y n +i(x ) = y 0 + f(t,y n (t )) dt 

xo 

produces a sequence of functions {y n (x)} that converges to this solution uniformly on I. 


Example 1: Consider the IVP 


y' = 3y 213 , 2 /( 2 ) = 0 

Then f(x,y) = 3 y 2 ^ 3 and = 2 y~ 1 ^ 3 , so f(x,y) is continuous when y — 0 but is not. 
Hence the hypothesis of Picard’s Theorem does not hold. Neither does the conclusion; the 
IVP lias two solutions, y 1 / 3 = x — 2 and y = 0. 


There are many ways to prove the existence of a solution to an ordinary differential 
équation. The simplest way is to ûnd one explicitly. This is a good approach for séparable 
or exact équations, or linear équations with constant coefficients. But unfortunately there 
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are many équations that cannot be solved by clcmentary methods, so attempting to prove 
the existence of a solution with this approach is not at ail practical. An alternative approach 
is to approximate a solution to an IVP by constructing a sequence of functions that converges 
to a solution. This is precisely the approach we will use for the proof of Picard’s Theorem. 

Before we discuss the idea bchind successive approximations, let’s first express a hrst- 
order IVP as an intégral équation. For the IVP y' = f(x,y), y(x o) = yo, suppose that / is 
continuons on sonie appropriate rectangle and that there is a solution y(x) that is continuons 
on some interval I. Then we may intégrale both sides of the DE to obtain intégral équation: 

X 

y{x) = y 0 + j f(t, y(t)) dt 

XQ 

Thus, under the assumptions of existence and continuity, the IVP is équivalent to the intégral 
équation. This fact seenis convenient and useful at first glance, but upon doser inspection 
we notice two problems: 


• The intégral équation is not wcll-defined unless we know that a solution exists. 1 

• The intégral équation is very hard to solve, except for very elementary IVPs. 


Suppose we dehne an operator T that maps a function y(x) to a function T[y](x), given 
by 

X 

T[y}(x) :=y 0 + j f(t,y(t)dt 

xo 

Then the intégral équation is simply y = T[y \, and any solution to the IVP must be a fixed 
point of T. 2 

To hnd fixed points, approximation methods are often useful. See Figure 1, below, for 
an illustration of the use of an approximation method to hnd a ûxed point of a function. 
To hnd a hxed point of the transformation T using Picard itération, we will start with the 
function yo(x) = y 0 and then iterate as follows: 


X 

y n +i(x) = y n (x) + J f(t,y n (t))dt 

XQ 

to produce the sequence of functions yo(x), yi(x), 1 / 2 (v), • • •• If this sequence converges, the 
limit function will be a hxed point of T. 

lr The terni ‘well-defined’ means essentially that ( t,y(t )) must be in the domain of f(x, y). For example, 
the function f(x) = is well-defined for ail x 7 ^ 1. But if x = 1, then f(x) is undefined. 

2 A fixed point of an operator or a transformation is an element in the domain that the operator or 
transformation maps to itself. In ternis of functions, a point x = a € dom g is a fixed point of g iff g (a) = a. 
Graphically, a fixed point is a point where the graph of y = g(x) intersects the straight line y = x. 
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Example 2: Consider the IVP 

y' = 2y, y{ o) = i 

x 

This IVP is équivalent to y = 1 + f 2ydt, so the Picard itérâtes are yo(x) = 1, 

o 


X 

Vi(x) = 1 + J 2y 0 (t) dt = l + 2x 

o 


y%(x) — 1 + J 2(1 + 2 1) dt — 1 + 2x + 


(2x) 

2 ! 


and so on. It can be shown by induction that the nth iterate is 

, . (2a;) 2 (2x) n " (2x 

Vn(x) = l + 2x+ { -^ + ...+ y -J- = J2 { 


2 ! 


which is the nth partial snm of the Maclaurin sériés for e 2x . Thus, as n —> oo, y n {x 


ni 

2x 


i —1 


l\ 


Jlx 


To carry ont a rigorous test for convergence (which is especially necessary when we don’t 
recognize the seqnence of Picard itérâtes), we need sonie idea of distance between functions. 
The distance measure used in the proof of the Picard Theorem is based on the norm of a 
function, as given in the following dehnition. 


Définition Let C[a,b\ dénoté the set of ail functions that are continuons on [a, b]. If y G 
C[a, b], then the norm of y is ||y|| := max \y(x)\. 

x£[a,b\ 

The norm of a function y(x) may be regarded as the distance between y[x) and y = 0, the 
function that is identically zéro. 3 With this in rnind we may define the distance between 
two functions, y,z G C[a, 6] to be the norm of y — z, or 


I \y-z 


max | y{x) — z(x) |. 

x£.[a,b] 


Using this measure of distance, we can define convergence of a sequence of functions to 
a limiting function. 

Définition A sequence {y n (x)} of functions in C[a,b\ converges uniformly to a function 
y(x) G C[a, b] iff 1 im \\y n - y\\ = 0. 


To illustrate uniform convergence, recall Example 2, where we found that y n (x) = Yi 

i=0 *' 

and y{x) = e 2x on the interval [0,1]. Then 


lim \\y n -y\\ = hni 


lim 

max 

y (2 xf e2 . 

= lim 

max 

^ (2 xf 


n—>oo 

xe[o,i] 

h 

oo 

xe[o,i] 

^ j ! 

i=n+1 L ' 



OO c\% 

= lim — = 0 

n—>oo L ' j ! 
i=n +1 L ' 


3 See Nagel, Saff, and Snider, page 837, for properties of this norm. 
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Therefore, y n —> y. One important detail to note in this example is that the uniform 
convergence of the sequence {y n {x)} to y(x) = e 2x on [a, b] occurs only when the interval is 
bonnded on the right. In other words, b mnst be hnite. If we were to take b to be infinité— 
for example, if the interval nnder considération were the whole real line—then the seqnence 
would not converge uniformly. The reason is that, if the value of x is nnbonnded, then for 

OO /y \ j 

any hnite n, the norm of X) y is inhnite. 

i=n-\-1 

Uniform convergence is particnlarly usefnl in that if a sequence of différentiable (and 
therefore continuons) funétions is uniformly convergent, then the function to which it con¬ 
verges is also continuons. 4 


Banach Fixed Point Theorem for Operators 

Let S dénoté the set of continuons functions on [a, b] that lie within a fixed distance a > 0 of 
a given function y t (x) G C[a, b], i.e. S = {y G C[a , b] : ||y — y t \\ < a}. Let G be an operator 
mapping S into S and suppose that G is a contraction on S, that is 

3k G R, 0 < k < 1 s. t. ||G[w] — G[z\ || < k\\w — z||Vu;, z G S. 

Then the operator G has a unique fixed point solution in S. Moreover, the sequence of 
successive approximations defined by y n+ \ := G[y n ], n = 0,1,2... converges uniformly to 
this fixed point, for any choice of starting function y 0 G S. 


To prove the theorem we ûrst show that the functions in the sequence y n + 1 = G[y n ] are 
wcll-dehned; that is, every function in the sequence {y n (x)} is in the set S. Next we show 
that this sequence converges uniformly to a function y^ G S. Finally, we show that this 
limit is a ûxed point of G , that is, y^ = Gly^]. 


Proof: 

Take any starting function y 0 G S. Since yo G Domc , then y\ = G[yo} is deûned. Since G 
maps S to itself, yi G S. By induction, y n G S and G[y n ] is well-deûned, for ail n > 0. 

Now rewrite y n = y 0 + {y± - y 0 ) + (y 2 - y i) + • • • + (y n - y n - 1 ), so that 


n— 1 

Vn{x) = y 0 (x) + - Vj( x )) (!) 

3=0 

We now show that the sequence {y n } converges uniformly to an element in S. We do this 
by applying the Weierstrass M-Test, an extension of the Comparison Test. 

4 See the text Introduction to Analysis by James R. Kirkwood, pages 206-212, for the définitions and 
proofs of some properties of uniform convergence. 
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Weierstrass M-Test Let {/„} be a sequence of functions defined on a set E. Suppose that 
for ail n G N, there exists M n G R such that 

\fn(x)\ < Mfix G E 

Then if J] M n converges, E fn must converge nniformly on E. 


What we need to do, then, is to find a bonnd M on the ternis (actually functions, of 
course) of the sériés (1). 


Claim: || y j+1 - y 3 ; || < k 3 \\y x - y 0 \\. 

The claim is evidently true for j = 0. Suppose that it is true for j = q, where q G N, q > 0. 

Then 

l|G[ Sa+2 ] - G[ Sa+1 ]|| = ||G[G[î,, +1 ]] - G[G[î,,]]|| < fc||G[ s , +1 ] - G[ s ,]|| < fc« +1 || Sl - So ||, 

proving the claim. 


Returning to the sériés (1), it is clear from the claim that max \y J+ \(x) — y 3 {'x) = \\y 3 +\ — 

x£[a,b] 

oo oo 

VjW < k 3 \\yi I/o||• Let M 3 := k 3 \\y x — 2/o||- Because E Mj = \\y 0 - y i|| E k 3 converges 5 , 

3 = 1 3=1 

the Weierstrass M-Test shows that {y n } converges uniformly to a continuons function, y^. 
Moreover, y^ G S because the assumption that Wy^ — y l \\ > a implies that \\y n — y t \\ > a 
for some n, contradicting the fact that y n G S. 

Since G is a contraction, we hâve that ||G[yoo] — G[t/ n ]|| < kWy^ — y n \\ for any n. But 
\\Voo ~ y n || —> 0 as n —> oo, so ||t/oc - y n+1 \\ —> 0 as n —> oo. Of course, G[y n ] = y n+1 . 
Thus, Jim ||G[t/oo] - G[y n \ || = ^lini^ ||G[yoo] - y n+1 \\ < ^lini^ fc||G[t/oo] - y n+ i\\ = 0. Finally, 

G[y oo] - y oc = (G{ yoo ] - y n+1 ) + (y n+1 - y^), so that 6 

||G[î/oo] VooW — ||G[|/oo] 2/n+l|| T llï/n+l 2/oo 1 1 • 

Since both ternis on the right side of the above équation approaeh zéro as n approaches oo, 
it follows that ||G[î/oo] — j/oo|| = 0, or G[yoo] = y^. Thus, t/oo is a ûxed point of G. 

Now suppose that z G S is any ûxed point of G , i.e. that z satisfies G[z\ = z. Then 
\\Voo ~ z\\ = ll G [l/oo] - G[z]\\ < kWyoo - z\\ < 11 y os - z\\, which is possible iff Wy^ - z\\ = 0. In 
other words, z = y^, so that y^ is the unique ûxed point of G. This complétés the proof of 
the Banach Fixed Point Theorem for Operators. 

Now we return to Picard’s Theorem. Suppose that u(x) is a solution to the IVP on 
[xo — h,x o + h] and recall that | f(x, y) < M on a rectangle Ri, as shown in Figure 2, below. 

5 It is a convergent géométrie sériés, because of the assumption that 0 < k < 1. 

6 By the Triangle Inequality for Norms. 
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Since u(x o) = yo, the graph of u(x) must lie in Ri for ail values of x close enough to xq. For 
such x , we hâve that \f(x,u(x))\ < M. Then \u\x)\ = \f(x,u(x))\ < M. So for x close to x$ 
the graph of u = T[u\ must lie within the shaded sector of Figure 2, where the parameters 
oq and h 2 are deûned geometrically. 

The graph of u(x) cannot escape from this rectangle on [rco — h 2 ,x o + h 2 \ since, if it 
did, \u'(x)\ = \f(x,u(x))\ > M at sonie point of the rectangle, which is clearly not possible. 
Thus, | u(x) — 2/0 1 < ot î, for ail x G [xq — h,xo + h]. So u(x) G S. 


Proof of Picard’s Theorem: 

To prove Picard’s Theorem we apply the Banach Fixed Point Theorem for Operators to the 
operator T. The unique ûxed point is the limit of the Picard Itérations given by 

y n+ i = T[y n ], y 0 (x) = y 0 . 


Recall that if y is a fixed point of T, then y(x) — yo+ f f(t,y(t )) dt, which is équivalent to 

xo 

the IVP. If such a function, y(x), exists, then it is the unique solution to the IVP. ' 

To apply the Banach Fixed Point Theorem for Operators, we must show that T will 
map a suitable set S to itself and that T is a contraction. This may not be true for ail real 
x. Also, our information pertains only to the particular intervals for x and y referred to the 
hypothèses of Picard’s Theorem. 

First we hnd an interval I — [x 0 — h,x 0 + h] and a G R, a > 0 such that T maps 
S = {g G C[I] : \\y— 2/0 1| < «} into itself and T is a contraction. Here, C[I) = C[xo — h,xo+h\, 
and we adopt the norm IMI/ = max \y\. Choose h\ and ai such that 

x£l 

Ri := {(x, y) : \x - x 0 \ < hi, \y — y 0 \ < ai} Ç R 

Because / and are continuons on the compact set Ri, it follows that both / and attain 
their supremum (and inhmum) on R ,\. 7 8 

It follows that there exist M > 0 and L > 0 such that 

V(x,y) G Ri, \f(x,y)\ < and 


7 Note that if there is a fixed point y(x), then y(x) = yo+ J f(t , y(t)) dt. By FTOC, y(x) is différentiable 

XQ 

and y'(x) = ^ (yo + f^ g f(t,y(t)) dt^J = f(x,y). Thus, y(x) satisfies the DE and solves the IVP. 

8 A compact set is usually defined to be a set with the property that if the set is covered by the union of 
ail members of a collection of open sets, then it is also covered by the union of a finite number of the open 
sets—any open cover has a finite subcover. An easier way to understand compactness is via the Heine-Borel 
Theorem and its converse which, together, state that a set in a Euclidean space is compact if and only if it 
is closed and bounded. For example, the set [0,1] is compact, since it is closed (i.e. contains its endpoints) 
and is clearly bounded. The set (0,1), on the other hand, is not compact. 


df 

dy 


< L. 
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Now let g be a continuons function on if = [x 0 — h\,xo + hi] satisfying \g(x) — yo\ < oy for 

X 

ail x G I. Then T[g\(x) = yo + J f(t,g(t )) dt so that, for ail x G I, 

XQ 


™W-ÿol = 

X 

j f(t,g(t))dt 

< 

X 

1 \ f(t,g(t))\dt 

< Ad 

X 

l ^ 


XO 


XQ 


XO 


M\x — x 0 


The situation is shown geometrically in Figure 2, below. 


Choose h such that 0 < h < minliq, ^}. Let a = oy, I = [x 0 — h,x 0 + h], and 
S = {g G C(J) : ||g — i/o||/ < a}. Then T maps S into S. Moreover, T[g](x) is clcarly a 
continuons function on I since it is différentiable (by the FTOC), and differentiability implies 
continuity. 


For any g G S, we hâve for any x G I, 

| T[g](x) - y 0 \ < M\x - x 0 \ < Mh < M (^j 

In other words, || T [g] — y 0 \\ < ai, so T [g] G S. 

Now we show that T is a contraction. Let u,v G S. On R\, 
there is a function z(t ) between u(t) and v{t) such that 


— oy. 


d£ 

dy 


< L, so by the MVT 


|T[«](*)-T[t;](x)| = 


< 


/ {f(t,u(t)) - f(t,v(t))}dt 

XQ 


j d / y (t,z(t)Mt)-v(t)\dt 



u(t)| dt 


< L\\u — v\\i\x — xo| < Lh\\u — v\\i, 


for ail x G /. Thus ||T[u](x) — T[n](x)|| < k\\u — n||, where k = Lh < 1, so T is a contraction 
on S. 


The Banach Fixed Point Theorem for Operators therefore implies that T lias a unique 
fixed point in S. It follows that the IVP y' = f(x,y),y(x o) = yo has a unique solution in S. 
Moreover, this solution is the uniform limit of the Picard itérâtes. 

Now we hâve found the unique solution to the IVP y' = f(x,y),y(x o) = yo in S, there 
is one important point that remains to be resolved. We must show that any solution to the 
IVP on I = [x 0 — h, x 0 + h] must lie in S. 

Suppose that u(x) is a solution to the IVP on [x 0 — h, x 0 + h]. Recall that | f(x, y)\ < M 
on the rectangle R\. Since u[x o) = yo, the graph of u(x) must lie in R\ for x close to xo- 
For such an x, we hâve that \f(x,u(x)\ < Ad, which implies that \u\x)\ = \f(x,u(x))\ < Ad. 
Therefore, for x close to Xo, the graph of u — T[u\ must lie within the shaded région of 
Figure 2. Moreover, the graph cannot escape from this région in [xq — h, Xq + h], since if it 
did, \u\x)\ = \f(x,u(x)\ > Ad at sonie point of the région, which is clearly impossible. Thus 
| u(x) — 2 / 0 1 < oy for ail x G [xq — h, Xq + h], which shows that u(x ) G S. 



